8.04 Quantum Physics 



Lecture IV 



Last time 

• Heisenberg uncertainty AxApx > | as diffraction phenomenon 

• Fourier decomposition 



dp(f){p)e'P''/^ (4-1) 
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V27Th 

m = -^m (4-2) 



Today 

• how to calculate (pik) 



• interpretation of iplx) and (pik): 

probability amplitude & probability density 



• measurement 



The larger (f){k), or (f){p), the more the wavefunction ip{x) resembles the plane wave, 
e^^^, that has definite momentum p = hk, and the more the particle described by the 
wavefunction iIj{x) is likely to be found to have momentum hk, if the momentum of the 
particle is measured. Conversely, if the particle's momentum is exactly p = po = hko, 
then the particle's wavefunction must be ip{x) = e*''"^, and the particle is equally likely 
to be found anywhere in space. Ax oo. To localize a particle in space, we need to 
add other Fourier components close to /cq. 




destructive interference 



Figure I: Addition of Fourier components close to ko yields a wavepacket localized in space. 
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To confine a wavefunction to a small region Ax of space, one needs many Fourier 

components, i.e., many plane waves of different momentum k G [ko, ko + Ak]. We will 
prove mathematically that, with suitably defined uncertainties Ax and Ak: 

follows strictly from the Fourier 

decomposition. Using Ap — UAk (4-3) 
we then arrive at 

Heisenberg uncertainty relation (4-4) 



Prom this viewpoint, the Heisenberg uncertainty relation follows from the decompo- 
sition of a wave into plane waves e*'^^, i.e., waves with definite wavevector /c, and from 
the relation between the wavevector k and the momentum p = hk. The smaller the 
region Ax, the more Fourier components p = hk are necessary to produce destructive 
interference everywhere outside Ax. 

Corollary: motion of particles. From the motion of particles and plane waves 
from Fig. I, it follows that if you change the relative phase of the different Fourier 
components, the constructive interference will occur somewhere else in space. This 
makes sense. We should be able to place a particle at a different location in space 
while maintaining the momentum distribution. It follows that if the relative phase 
between the plane waves changes continuously, the location of constructive interfer- 
ence (i.e., the location of the particle) will move in space. In terms of wave mechanics, 
the motion of the particle is simply due to a change of phase between the Fourier 
components (i.e., the plane waves). Then, to reproduce CM, the phase of the Fourier 
components must rotate in time at a frequency that depends on the momentum p. 

Corollary: time-energy uncertainty. If instead of Fourier transforming the posi- 
tion coordinate at fixed time \E'(a;, to) ^{k, to) we fix the position X — Xq and study 
the time variation of the wavefunction ^(xo,^), then we can Fourier decompose the 
wavefunction into frequency components. 



AxAk > - 



AxAt) > — 
^ - 2 



dujf{uj)e 



-iu)t 



(4-5) 



Convention. Positive frequency u corresponds to negative phase evolution. Applying 
the same mathematical and logical arguments as before, we arrive at: 



AuuAt > 



1 



AEAt > 



h 



time-frequency uncertainty 

energy-time Heisenberg un- 
certainty relation {E — fiuj) 



(4-6) 
(4-7) 
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V(Xo,t) 




Figure II: Wave function vs. time. 



Probability amplitude and probability density 

For light or other waves, the energy per unit volume (or per unit length) is propor- 
tional to the square of the electric field. Since the number of photons per unit volume 
is proportional to |Ep, we postulate in analogy: 



|-?/'(x)p is called the probability density (probability per unit length). The waveform 



Figure III: The magnitude squared |^/^(x)p of the wavefunction is the probabihty density 
to find the particle in a region of space. The probability to find the particle in the spatial 
interval {x,x + dx] is given by \ilj{x)\'^dx. 

iIj{x) is also called the probability amplitude (more exactly: probability density ampli- 
tude). In contrast to EM fields, \E' is truly a complex quantity. The requirement that 
the particle be found somewhere in space leads to the normalization condition: 




(4-8) 




IV(x)P 



X 




(4-9) 



In the homework, you will prove Parseval's theorem: 



Proof. If <f){k), (f){p), is the Fourier transform of a wavefunction iplx), i.e., if 
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then, 

/oo poo poo 

dx\i){x)\'^ = / rfA;|0(A;)|2= / c/p|0(p)p. □ 
-oo J —oo J ~oo 

It follows that if the wavefunction |'?/'(a;)p is normalized, so is We already 

argued that if (j){p) is peaked around some value po? the motion of the particle will 
be similar to that of a classical particle with momentum po? (plane wave e*^"^'/^). 
Taking into account Parseval's theorem, it is reasonable to interpret (f){p) as the 
probabilty amplitude for momentum, i.e., 

/ Probability to find a particle momentum \ _ . , . ,2 , / . ,„n 

y between p and p + dp J \ \ )\ ■ 

Similarly, is the probability density for the wavevector k. Note that Aa; (or 




p p+dp 



Figure IV: The probability density in momentum space is given by 



Ap) is not the uncertainty of the measurement device, but associated with the particle 
itself. If our measurement apparatus can resolve with resolution Axapp <^ Ax, and 
we repeat the experiment with an identically prepared particle many times, we 
will observe a histogram. In the limit of a very large number of measurements, the 
histogram reproduces the probability density |-?/^(x)p or |0(j9)p. 




\ >i Pi 



Figure V: Probability density recon- Figure VI: . . . and in momentum 
struction in position space. space. 

Note. After measuring a particular value x with apparatus uncertainty Axapp, the 
particle is no longer described by the original wavefunction ^/'(a;), but by a new wave- 
function ipi^x) that is consistent with the outcome of the measurement result ("collapse 
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of the wavef unction"). In particular, if Axapp -C Ax, the spread of the new wave- 
function ip{x) in momentum wiU be much larger than before, consistent with the 
Heisenberg uncertainty {Ax = Ax^ = Axapp). 

AxAp = AxappAp > |, (4-11) 



or 



A^ > 



n 



2Ax 



> 



app 



2Ax 



2' 
= Ap. 



(4-12) 



If you now choose to measure momentum with resolution Apapp ^ Ap, the uncer- 
tainty in position will again increase and so on. 



■ip{x) 



1 



dk(j){k)e 



ikx 



dp4>{p)e 



ipx/h 



(4-13) 



How to determine the momentum distribution given 
the wavefunction iIj{x) in position? 

The expansion coefficients (f){k) are given by the inverse of the Fourier decomposition: 



m 



'2tx J 



dxip{x)e 



-ikx 



Proof. 



dk(f){k)e 



ikx 



1 

2^ 
1 

2^ 



ikx^ikx 



dk I dx'ip{x')e 
dx'i){x') I c/A;e^^(^-^') 



□ 



What is the value of dke^^^^~^'^l Qualitatively, if x 7^ x', the integrand oscillates 
in the complex plane many times a A; — > ±cxd, so the integral is zero. If x = x', 
the integral is dk ■ 1 and diverges. So the "function" I{y) = dke^^^ looks 
something like: 



I(0)=oo 


i(y) 


I(yy^O)=0 


I(y7^0)=0 




Figure VII: Sketch of the "function- 
represented by the integral I{y). 



Figure VIII: Convolution of I{y) 
and a Gaussian function to deter- 
mine "the area under I{y)" ■ 
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How bad is the divergence? Let us calculate the area under the curve a real and 
positive. 

/oo poo poo 

dyl{y)e-''y' = / dy dke'^^ e'^'y'' 
OO J —OO J —oo 

/oo roo 
dk / dye-'^y^^'^y (4-14) 
-oo J — oo 

To calculate the integral, we note without proof that: 

j dye~°'^y~^^'^ — for any complex a, {3 with Re(Q;) > 0. (4-15) 

To bring the above integral, Eq. (4-14), into the desired form, we expand the expo- 
nent: 

-ay +zky^-a^y --y+[^) j-^ 



/oo /*00 /*CX) 

dyl{y)e-'^y" = dk dye""!^-^) 
■oo J —oo J —oo 



^ 1 oo ., 



dke 4a 



Vtt • 4q; 

= 27r (4-17) 
independent of the value a! Since our probing function e~"^ > 1, we con- 

a—^0 

elude that the area under our function I{y) is finite and equal to 27r. We define a 
"generalized function" (mathematically, a distribution) by 

5{x) = l' (4-18) 
\oo for a; = 0. 

with the property, 

/oo 

dxS{x) = 1. (4-19) 
-oo 

This is the Dirac delta function. We can think of it as the limiting case of a 
function of finite width (e.g., a Gaussian or square function) that is made narrower 

and narrower, while keeping the area under it constant. We have dke^^^ — 27rS{y). 
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Figure IX: The delta function Figure X: Dirac delta function, 

can also be expressed as 5{x) = 



,2 



Properties of the delta function 

What is dxf{x)6{x — XQ)'! For a function f{x) that is sufficiently smooth (regular) 
at X = and using 6{x — xq) = for x 7^ xq, we have 



dxf{x)6{x — Xq) 



XQ+t 



dxf{x)6{x — Xq) 



XQ-e 



XQ+e 



dx6{x — Xq) 



xo-e 
00 



f{xo) / dx6{x) = f{xo). 



(4-20) 



Therefore, we have: 



dxf{x)5{x - Xq) = f{Xo) 



(4-21) 



Convolution of a function f{x) with 5{x — Xq) "projects out" the value of the function 




Figure XI: Convolving a sufficiently smooth function with a Dirac delta function projects 
out the function value at one point. 
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e(x) 



Figure XII: The 5-function is also the derivative of the Heavyside step function. 



at Xq. Without proof, we note that S{x) = ^B(x). Derivative of the deha function 
can be defined. Integration by parts yields 

J dxf{x)5'{x - xo) = -/'(x - xo). (4-22) 

Convolution with 5' projects out the negative derivative at xq. 
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